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Abstract 

We present the next-to-next-to-leading order (NNLO) corrections to the 
total cross section for (pseudo-) scalar Higgs boson production using an alter- 
native method than those used in previous calculations. All QCD partonic 
subprocesses have been included and the computation is carried out in the 
effective Lagrangian approach which emerges from the standard model by 
taking the limit — oo where mt denotes the mass of the top quark. Our 
results agree with those published earlier in the literature. We estimate the 
theoretical uncertainties by comparing the ii'-factors and the variation with 
respect to the mass factorization/renormalization scales with the results ob- 
tained by lower order calculations. We also investigate the dependence of 
the cross section on several parton density sets provided by different groups. 

-'^partially supported by the National Science Foundation grant PHY-0098527. 



Further we study which part of the coefficient functions dominates the cross 
section. This is of interest for the resummation of large corrections which 
occur near the boundary of phase space. It turns out that depending on the 
definition of the total cross section the latter is dominated by the the soft- 
plus- virtual gluon corrections represented by S{l—x) and (ln*(l— 
terms. 

PACS numbers: 12.38.-t, 12.38.Bx, 13.85.-t, 14.80.Gt. 
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1 Introduction 



The Higgs boson, which is the corner stone of the standard model, is the only 
particle which has not been discovered yet. Its discovery or its absence will 
shed light on the mechanism how particles acquire mass as well as answer 
questions about supersymmetric extensions of the standard model or about 
compositeness of the existing particles and the Higgs boson. The LEP exper- 
iments [1] give a lower mass limit of about mn ~ 114 GeV/c^ and fits to the 
data using precision calculations in the electro-weak sector of the standard 
model indicate an upper limit itlh < 200 GeV/c^ with 95 % confidence level. 
After the end of the LEP program the search for the Higgs will be continued 
at hadron colliders in particular at the TEVATRON and the LHC. 

In this paper we concentrate on Higgs production where the lowest order 
reaction proceeds via the gluon-gluon fusion mechanism. In the standard 
model the Higgs boson couples to the gluons via heavy quark loops. Since 
the coupling of the scalar Higgs boson H to a fermion loop is proportional to 
the mass of the fermion (for a review see [2]), the top-quark loop is the most 
important. The latter contribution is also dominant for the pseudo-scalar 
Higgs boson A provided tan/3 is small where (3 denotes the mixing angle 
in the Two-Higgs-Doublet model (2HDM). In lowest order (LO) the gluon- 
gluon fusion process g + g ^ B with B = H, A, represented by the top-quark 
triangle graph, was computed in [3]. The two-to-two body tree graphs, given 
by gluon bremsstrahlung g+g — > g+B, g+q{q) — > q{q)+B and q+q — > g+B 
were computed for B = H in [4] and for B = A in [5]. From these reactions 
one can derive the transverse momentum (px) and rapidity (y) distributions 
of the (pseudo-) scalar Higgs boson. The total integrated cross section, which 
also involves the computation of the QCD corrections to the top-quark loop, 
has been calculated in [5]. This calculation is rather cumbersome since it 
involves the computation of two-loop triangular graphs with massive quarks. 
Furthermore also the two-to-three parton reactions have been computed in [6] 
using helicity methods. From the experience gained from the next-to-leading 
(NLO) corrections presented in [5] it is clear that it will be very difficult to 
obtain the exact NLO corrections to one-particle inclusive distributions as 
well as the NNLO corrections to the total cross section. 

Fortunately one can simplify the calculations if one takes the large top- 
quark mass limit mj — > cxo. In this case the Feynman graphs are obtained 
from an effective Lagrangian describing the direct coupling of the (pseudo-) 
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scalar Higgs boson to the gluons. The LO and NLO contributions to the total 
cross section in this approximation were computed in [7] . A thorough analysis 
[5], [8] reveals that the error introduced by taking the — > oo limit is less 
than about 5% provided uih < 2 rrit. The two-to-three body processes were 
computed with the effective Lagrangian approach for the scalar and pseudo- 
scalar Higgs bosons in [9] and [10] respectively using hclicity methods. The 
one-loop corrections to the two-to-two body reactions above were computed 
for the scalar Higgs boson in [11] and the pseudo-scalar Higgs boson in [12]. 
These matrix elements were used to compute the transverse momentum and 
rapidity distributions of the scalar Higgs boson up to NLO in [13], [14] and the 
pseudo-scalar Higgs boson in [12]. The effective Lagrangian method was also 
applied to obtain the NNLO total cross section for scalar Higgs production by 
the calculation of the two-loop corrections to the Higgs-gluon-gluon vertex in 
[15], the soft-plus- virtual gluon corrections in [16], [17] and the computation 
of the two-to-three body processes in [18], [19]. These calculations were 
repeated for pseudo-scalar Higgs production in [20], [21]. 

In this paper we recalculate the NNLO corrections to the total cross sec- 
tion for (pseudo-) scalar Higgs production which have been computed recently 
in [18]- [21] and we found complete agreement with their analytic results when 
the number of colours is set N = 3. However our method differs from the one 
presented in [18], [20] and the approach followed in [19], [21]. The authors 
in the latter references compute the total cross section using the Cutkosky 
[22] technique where one- and two-loop Feynman integrals are cut in certain 
ways. These Feynman integrals can be computed using various techniques 
(for more details see the references in [19]). Furthermore this method been 
also applied to compute rapidity distributions in NLO which is presented 
in [23]. The authors in [18] choose a more conventional method which was 
already used in [24], [25], [26] to compute the coefficient functions for the 
Drell-Yan process. However instead of an exact computation of the 2 — > 2 
and 2 — » 3 body phase space integrals they expand them around x = 1. Here 
X = / s where m and ^/s denote the (pseudo-) scalar Higgs mass and the 
partonic centre-of-mass energy respectively. Since the coefficient functions 
can be expressed into a finite number of Polylogarithms of the types Li„(a;), 
^n,p{x) (see e.g. [27]) and logarithms of the types In'x ln-'(l — x), which 
are all multiplied by polynomials in x, one can expand them in the limit 
X ^ 1 and match them with the expressions coming from the phase space 
integrals. In this way the coefficients of the above functions are determined. 
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However as is shown in [26] this procedure is not necessary since one can 
obtain the Polylogarithms in a direct way by expanding the hypergeometric 
functions [28] in e — n — A where n is the number of dimensions characteristic 
of n-dimensional regularization. The point is that there is no essential dif- 
ference between Higgs production in the effective Lagrangian approach and 
the Drell-Yan production mechanisms except that the matrix elements in the 
former case contain higher powers of the invariants in the numerators. In 
[26] a program was made using the algebraic manipulation program FORM 
[29] to evaluate the phase space integrals analytically. For this calculation 
we have extended the program so that it can accommodate integrals having 
higher powers of invariants in the numerator. Our approach can be also used 
for differential distributions in particular for jet production (see e.g. [30]). 
Here one cannot avoid multi-particle phase space integrals. Of course two- 
to-four body processes are even more cumbersome to deal with but one can 
at least try to compute the pole terms {1/e)^, which arise from infrared and 
collinear divergences, analytically or numerically (for numerical methods see 
[31]). Another difference with the previous results is that we present the 
radiative parts of the coefficient functions for general colour factors of the 
local gauge group SU(N). The latter is important if one wants to resum the 
dominant contributions, (see [8], [32]). 

Our paper will be organized as follows. In Section 2 we give definitions 
and discuss the regularization method to compute the partonic cross sections 
in the effective Lagrangian approach. In Section 3 we present the calculation 
of the coefficient functions. In particular we discuss the various frames in 
which the three-particle phase space integrals are computed. In Section 4 we 
present the total cross sections for scalar Higgs production at the LHC and 
the TEVATRON. The long expressions for the coefficient functions expressed 
in the various colour factors are presented in Appendix A (scalar Higgs) and 
Appendix B (pseudo-scalar Higgs). 
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2 Application of the effective Lagrangian ap- 
proach to Higgs production 

In the large top-quark mass limit the Feynman rules (see e.g. [9]) for scalar 
Higgs production (H) can be derived from the following effective Lagrangian 
density 



£H ^ = Gh <^^{x) 0{x) with 0{x) = 
whereas pseudo-scalar Higgs (A) production is obtained from 



-Gl^{x)G'^'^'^{x), (2.1) 



GaOi(x) + Ga02(x] 



with 



Oi(x) = -^e^.A.G'rG'f(x), 



^ i=i 



(2.2) 



where $^(x) and $^(x) represent the scalar and pseudo-scalar fields respec- 
tively and Uf denotes the number of light flavours. Furthermore the gluon 
field strength is given by G^" and the quark field is denoted by g^. The fac- 
tors multiplying the operators are chosen in such a way that the vertices are 
normalised to the effective coupling constants Gh_, Ga and Ga- The latter 
are determined by the top-quark triangular loop graph, including all QCD 
corrections, taken in the limit — > oo which describes the decay process 
B ^ g + g with B = H, A namely 



Gb = -2^/^a,(/x^)G]/VBFB(rB)CB(^a,(/x^),^ 



and as{n1) defined by 



2n «s(/^r) 



Ga : 



(2.3) 



(2.4) 
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where as(yU^) is the running couphng constant and fir denotes the renormal- 
ization scale. Further Gp represents the Fermi constant and the functions 
Fb are given by 

Fnir) = 1 + (1 - r) /(r) , FA(r) = /(r) cot /? , 
4 m? 



fir) 



arcsm 



for r > 1 , 



/w = -J 



In 



1 + ^ 



TT 2 



for T < 1 , 



(2.5) 



where cot f3 denotes the mixing angle in the Two-Higgs-Doublet Model. Fur- 
ther m and rrit denote the masses of the (pseudo-) scalar Higgs boson and 
the top quark respectively. In the large mrlimit we have 



lim Fnir) = — , 

r— »oo o T 



lim Fa(t) = - cot/?. 

r— »oo T 



(2.6) 



The coefficient functions Cb originate from the corrections to the top-quark 
triangular graph provided one takes the limit — > oo. We have presented 
the Born level couphngs Gb in Eq.(2.3) for general nit for on-shell gluons 
only in order to keep some part of the top-quark mass dependence. This is 
an approximation because the gluons which couple to the H and A bosons 
via the top-quark loop in partonic subprocesses are very often virtual. The 
virtual-gluon momentum dependence is neither described by Fb{t) nor by Cb- 
However for total cross sections the main contribution comes from the region 
where the gluons are almost on-shell so that this approximation is better 
than it is for differential cross sections with large transverse momentum. 
The coefficient functions are computed up to order in [8], [33] for the H 
and in [34] for the A. They read as follows 



V rntj 



hC A — 3 Cl 



2 ^ 



CATf + {lC^ 
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mj \ 9 mi 



(2.7) 



C,^{as{i4).^=l. (2.8) 

where af^ is presented in the five-flavour number scheme. Notice that the 
coefficient function in Eq. (2.7) is derived for general colour factors of the 
group SU(N) from Eq. (6) in [33]. These factors are given by 

N'^ - I 1 
Ca^N , Cf = ^^ , Tf = -. (2.9) 

Notice that T/ is also incorporated into Gb in Eq. (2.3) where it is set to the 
value Tf = 1/2. 

Using the effective Lagrangian approach we will calculate the total cross 
section of the reaction 

/7i(Pi) + H2iP2) ^ B{-p,) +' X' , (2.10) 

where Hi and H2 denote the incoming hadrons and X represents an inclusive 
hadronic state. The total cross section is given by 

° ^1 a,b=q,q,g ''^ ''^/^^ 

(X Tf^ \ 
,— , B = H,A, 
Xi X2 jJ.^ ) 

With x = ^ , 5 = (Pi + P2)' , vl-m\ (2.11) 

where the factor l/^N"^ — 1) originates from the colour average in the case of 
the local gauge group SU(N). Further we have assumed that the (pseudo-) 
scalar Higgs boson is mainly produced on-shell i.e. p\ ~ w?. The parton 
densities denoted by fa{y^fJ^) (a, & = q,(l,g) depend on the mass factor- 
ization/renormalization scale fi. The same scales also enter the coefficient 
functions Aab,B which are derived from the partonic cross sections 

^' "si (2^ ^^^'^^ ~ "^^^ Tab,B{p5,Pl,P2) , (2.12) 



where the incoming parton momenta, denoted by pi and p2, are related to 
the hadron momenta by 



m 



s = xiX2S, z~ — . (2.13) 

s 



S={pi+P2f , =^ 

The amphtude Tab,B can be written as 

TabAP^,Pi,P2) = Glld'ye'P-y{a,b\0(y)Oma,b), (2.14) 

TabAP5:PuP2) = Jd^ye'P^-y{a,b\(GA0,{y) + GA02{y)) 

X (Ga Oi(0) + Ga 02(0)) |a, b) . (2.15) 

The expressions above for the amphtude T^b are similar to those given for 
the Drell-Yan process except that the conserved electroweak currents are 
replaced by the operators O and Oi,02. The latter are not conserved so 
that they acquire additional renormalization constants defined by 

0{y) = Zod{y) , 0,{y) = Z,,dj{y). (2.16) 

where the hat indicates that the quantities under consideration are unrenor- 
mahzed. Insertion of the above equations into Eqs. (2. 14), (2. 15) leads to the 
renormalized expressions 

TabAP5:Pi:P2) = G^Zl [ ^ e^^^'M^, &|0(2/) 0(0) |a, 6) , (2.17) 



Tab,AiP5,Pi,P2) = J d^ye'P'-y 



^A ^11 + ^A + 2 Ga Ga Zii Z21 



X (a, 6|Oi(y) Oi(0) |a, 6) + Zl^ + Gl Zl^ + 2 Ga Ga ^12 ^22 



X (a, b\d2{y) 02(0) |a, 6) + G^ Zn Z12 + Gl Z22 ^21 + Gj, Ga 



X (Zn ^22 + ^12 ^21) [ (a, 6|Oi(y) 62(0) + ^(y) Oi(0)|a, 6) 



(2.18) 



The operator renormalization constants depend on the regularization scheme 
in particular on the prescription for the 75-matrix and the Levi-Civita tensor 
in Eq. (2.2). The computation of Tab,H will be carried out by choosing n- 
dimensional regularization where in the case of Tab,A we adopt the HVBM 
prescription [35], [36] for the 75-matrix. For this choice the contraction of 
the Levi-Civita tensors proceeds as 



§U2 §X2 ^CJ2 

"/il "/ii "/ii 

Kl ^il 

K ^xl K 

^U2 (5^2 



(2.19) 



Ul Ul (Tl 

where all Lorentz indices are taken to be n-dimensional. To facilitate the 
calculation one can replace 7^75 in Eq. (2.2) by (see [37]) 



(2.20) 



which is equivalent to the HVBM scheme. Choosing the MS subtraction 
scheme the renormalization constant corresponding to the operator O be- 
comes [38] 



£ 



4 2 

e 



+ 



(2.21) 



where Ss denotes the spherical factor characteristic of n-dimensional regular- 
ization. It is defined by 



expi- 



7e — ln47r 



(2.22) 



The lowest order coefficients and (3i originate from the beta-function given 

by 



34 



20 



rifTf, p,^ Ci-4nfTfCF--nfTfCA. 



(2.23) 
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The operator renormalization constants corresponding to Oi and O2 are com- 
puted in [39] and they read 



Zll = ^a. = 1 + dsil^r) Se-pO + O^l/^r) 

e 



+ 



(2.24) 



e 

where Z^^ denotes the couphng constant renormahzation factor defined by 

ds = Za, as{i4) ■ (2.25) 
The remaining constants are 

Z21 = 0, (2.26) 



Z12 = Gsinl) Se-l-QCp 

e L 

^22 — ^MS ^5 ' 



(2.27) 

(2.28) 



where Z^g and Z| are the constants characteristic of the HVBM scheme. 



They are given by 



1 + as{iJil 



-4Cf 



2 / 2\ 



107 



31 



22Cj,-—CACF + —njT^C,^ 



9 



9 



(2.29) 



The latter renormahzation constant is determined in such a way that the 
Adler-Bell-Jackiw anomaly [40] 



(2.30) 



is preserved in all orders in perturbation theory according to the Adler- 
Bardeen theorem [41]. 
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3 Computation of the partonic cross section 
up to NNLO 

In this section we will give a short outline of the computation of the partonic 
cross sections aab,B in Eq. (2.12) up to next-to- next-leading order (NNLO). 
Insertion of a complete set of intermediate states and using translational 
invariance in Eq. (2.13) we obtain for scalar Higgs production 

° •' m=3 1=3,17^5 •' 

m 

J=l 

Mah^X H = (Pl,P2|O(0)|X,p5) with |X,P5) = |P3,P4,P6 " " ' Pm-/Pf,) , 

(3.1) 

where Kah represents the spin and colour average over the initial states. Fur- 
ther one can write down a similar expression for the pseudo-scalar Higgs (see 
Eq. (2.15)). Up to NNLO wc have to compute the following subprocesses. 
On the Born level we have the reaction 

g + g^B. (3.2) 

In NLO we have in addition to the one-loop virtual corrections to the above 
reaction the following two-to-two body processes 

g + g^B + g , g + q{q) ^ B + q{q) , q + q^B + g. (3.3) 

In NNLO we receive contributions from the two-loop virtual corrections to 
the Born process in Eq. (3.2) and the one-loop corrections to the reactions 
in Eq. (3.3). To these contribution one has to add the results obtained from 
the following two-to-three body reactions 

g + g^B + g + g , g + g^B + q^ + q^, (3.4) 
g + q{q) ^B + q{q)+g, (3.5) 
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q + q^B + g + g , q + q^B + qi + qi, (3.6) 

gi + 52 ^ B + gi + g2 , gi + ^2 ^ B + + ^2 , (3.7) 

q + q^B + q + q. (3.8) 

In the case of pseudo-scalar Higgs production one also has to add the con- 
tributions due to interference terms coming from the operators Oi and O2 in 
Eq. (2.18). Up to order we have to compute the following expression for 
the reactions in (3.4) 

^ •' LI m=3 i=3 •' K"^^) 



X 



{{Pl,P2\Ol\P5) {P5\02\Pl,P2) + ipi,P2\02\P5) {P5\Ol\Pl,P2 



. (3.9) 



The Feynman graphs for {pi,P2\Oi\p5) and {pi,P2\02\p5) can be found in 
Fig. la and Fig. 2a of [20] respectively. From the above equation one infers 
that the interference term contributes as a delta-function S{1 — z) to o"c,c,,a- 
Furthermore since is proportional to Ga (see Eq. (2.3)) one can extract 
the latter constant as an overall factor from the above equation. For the 
reactions in Eq. (3.5) we have to compute 



ab 



dy,5+(pl-m') 



' 4 m 

m=3 «=3 ^ C^^. 



d^-pi 



x5("nEp.)l(Pl,P2|Ol|^,P5)r| + |^12 / (2^'^^^^3) 
X5(")(pi-Fp2+P3+P5) ((Pl,P2|C>i|p5,P3) (P5,P3|02|P1,P2) 



+ (P1,P2\02\P5,P3) {p5,P3\Ol\Pl,P2 



(3.10) 
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The Fcynman graphs for {pi, P2\0i\pz, Pb) and (pi,P2|02|p3,P5) can be found 
in Fig. lb and Fig. 2b of [20] respectively. In the HVBM-scheme the latter 
matrix element is proportional to £ = n — 4. However this contribution does 
not vanish in the limit £ — > because of the single pole term present in Z12 
in Eq. (2.27). The same expression as in Eq. (3.10) also exists for (yqq^A 
originating from the first reaction in Eq. (3.6). In this case the Fcynman 
graphs corresponding to (pi,P2|Oi|]93,]95) and (pi,P2|02|p3,]?5) are shown in 
Fig. Id and Fig. 2d of [20] respectively. Notice that these contributions also 
survive in the calculation of differential distributions, see [12]. 

The computation of the phase space integrals proceeds as follows. First we 
take over the expressions from the one- and two-loop corrections to the Born 
reaction in Eq.(3.2) which are computed in [7] (one-loop) and [15] (two-loop). 
Unfortunately the two-loop result is not presented for arbitrary colour factors 
so that we cannot distinguish between n/ Ca Tf and UfCpTf. Hence we have 
shuffled all contributions proportional to Uf into the term proportional to 
UfCATf. To compute the two-to-two body processes including the virtual 
corrections we have chosen the centre-of-mass frame of the incoming partons 
given by 

Pi = i/^(l,0,---,0,l). 




P2 = -VA2(l,0,---,0,-l), 



2 



Pi2-m 

P3 — , (1,0, • • •, — sm^, — cos^J , 




-P5 = ;7-7== iPi2 + 0, • • • , (P12 - m^) sin 6, (P12 - m^) cos 6) , 

(3.11) 

where we have introduced the invariants 

P^, = iP^+PJy , Pl2 = S. (3.12) 

In this frame the two-to-two body phase space integral becomes 

^/-/2 93-n -1 /p \ l-"/2 / p 2\"-3 

_ ry2 ^*-^H ^ ^ I -'12 \ I -'12 — m \ 

Ca&-»c H — -"-a&^O 



P12 (47r)"/2 r(n/2 - 1) \ n'^ J \ fi 



2 
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Jo 



X 



d6 sin 



n—3 



e\M, 



2 



(3.13) 



ab^c B 



Here fi indicates the dimension of the strong couphng constant g ^ g /x^^"")/^ 
which is characteristic of n-dimensional regularization. In order to make all 
ratios dimensionless we have ah^cady included the factor in the phase 

space integrals. Notice that in principle the scale fi has nothing to do with the 
factorization or rcnormalization scale. However for convenience one puts the 
latter scales equal to ii. The evaluation of the integral in Eq. (3.13) is rather 
easy even when \Mab^c b| contains virtual contributions. Since the integrals 
can be evaluated analytically we have made a routine using the algebraic 
manipulation program FORM which provides us with the results. A similar 
routine is made for the two-to-thrcc body phase space integrals but here the 
computation is not so easy unless one chooses a suitable frame. Since we 
integrate over the total phase space the integrals are Lorentz invariant and 
therefore frame independent. The matrix elements of the partonic reactions 
can be partial fractioned in such a way that maximally two factors Pij in Eq. 
(3.12) depend on the polar angle 6i and the azimuthal angle 62- Furthermore 
one factor only depends on 9i whereas the other one depends both on 9i and 
62- Therefore the following combinations show up in the matrix elements 



For the first combination it is easy to perform the angular integration since 
all momenta represent massless particles and one obtains a hypergeometric 
function (see e.g. Eq. (4.19) in [14]). The angular integral of the second 
combination is more difficult to compute because one particle is massive and 
the result is an one dimensional integral over a hypergeometric function which 
however can be expanded around e. Examples of these types of integrals can 
be found in Appendix C of [42]. The last combination is very difficult to 
compute in n dimensions because both factors contain the massive particle 
indicated by ^5. Therefore one has to avoid this combination at any cost. 
This is possible if one chooses the following three frames. The first one is the 
centre-of-mass frame of the incoming partons. The momenta are 



yk pi pk pi pk pi 
ij mn ' ij m5 ' i5 m5 



4>A;>-2,4>/>-2, 



222 2 n 

Pi =Pj=Pm=Pn^^ 



2 2 

P5 = m . 



(3.14) 



Pi 
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P2 = -VA2(1,0,---, 0,0,-1), 



-P3 = 



-Pi 



(1, 0, • • • , 0, sin^i, cos^i) , 

(1, 0, • • • , sm%l) sin 621 cosip sin 9i + sinip cos ^2 cos^i. 



• 2 "0 
2 



-P12 — P45 

cos'^ cos^i — sin'^ cos ^2 sin^i) , 

P12 -P34 
{P12 - P35) (P12 - P45) ' 



(3.15) 



7^ 72 1 



^ Pu 

O r,lD / A^\n T^/„ 0\ \ ,.2 



l-n/2 



2Pi2 (47r)" r(n - 3) V /^^ 

P35Pi5-Pl2m' 



Pl2 



dP 



35 



Pl2+m2-P35 /D D D ^2\«/2-2 



X / (iP45 
'Pi2 mVPss 



X 



'Pl2+m2-P35-P45' 



n/2-2 



Jo 

X d^2 sin'^-' 02 |M„6^ed bI' ■ (3.16) 
Jo 

For the next frame we choose the centre-of-mass frame of the two outgoing 
particles indicated by the momenta and p^. 

pi = wi(l,0,---, 0,0,1), 

P2 = (a;2,0, • • • ,0, IpsI sinV', IpsI cosV' - <^i) , 
1 



-P3 



' P34 (1, 0, • • • , sin 9i sin 62, sin 61 cos 82, cos ^1) , 



1 



— P4 = - y P34 (1, 0, • • • , — sin^i sin^2, — sin^i cos^2, — cos^i) , 
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-P5 



{uj5,0, • • •, 0, |p5| sin^, |p5| cosi/j) , 



Pi2 + Pir, -m" Pu + - 

- , (^2 = ■ 



UJ5 = 



2VPu 

Pl5 + P25 



34 



COS'0 = 



2^ ' 

P34 - (Pi5 - m2) - P12 (P25 + m2) 



(P12 + Pi5 - m2) J(Pi5 + P25)2 - 4m2 P, 



(3.17) 



34 



rmHPi2 my{P2,-m^) ^ ^ f (p^^ _ ^2-) ^p^^ _ ^2) _ p^^ ^2 n/2-2 
Jn 



w?+Pi2 m?/{P25-m^) 
X / dPi5 



X r (102 8111^-^62 |M< 



ab-^cd B 



(3.18) 



For the last frame we choose the centre-of-mass frame of one of the outgoing 
partons and the (pseudo-)scalar Higgs boson indicated by the momenta 
and respectively 

pi = u;i(l,0,---, 0,0,1), 

P2 = (a;2, 0, 0,ci;3 sin Go's COS co-i) , 

-P3 = a;3(l,0, •••,0,sinV',cosV'), 

— P4 = a;4 (1, 0, ■ • • , sin6'i sin6'2, sin 6'i cos6'2, cos6'i) 

—p5 — {005, 0, • • • , —004^ sin 9i sin 02, —uja sin cos ^2, —(^4 cos ^1) , 
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Pi 
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2a/P 
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2^/P 
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45 



CU4 — 



45 



2^/P 



45 



<^5 = 



P45 + 
2a/P 



45 



COS'0 = 



Pi 2 P23 ~ P4R P 



45 -ri3 



(P12 + P13) (Pl3 + P23) ' 



(3.19) 



, ,p ,'Pl3P23\"^' V^12 + ^13 + P23-m2V ' 



X 



P12 + Pl3 + P23 \ 



7o 7o 



x|M„5-.cdBr. (3.20) 

The integration over the angular independent invariants Pij can be performed 
by rescahng them with respect to P12 (see Appendix E in [25]). The integrals 
are such that they can be performed in an algebraic way by a program based 
on FORM [29] which has been also used to compute the NNLO coefficient 
functions of the Drell-Yan process in [24], [25], [26] (for more details see [43]). 
The partonic cross sections Eq. (2.12) can be written as follows 

a„.,Hl^^,^J - 8(Ar2-i)TT^^o^'^^'" l^^vj 

/ nGl l + £ ( rn?\ 



m 



1^' 



e/2 



4°^-2/3o)+^Sb + ^^Sb 
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^ (2 pS> » PS' + » Pi? - 10 A Pif 



(8) w 



(1) 



+2 ^Sb ® n^? + ^b) - 6 /3o <;,B + 2 Pi? ® 



(1) 

Sfl>B 



with — —4 /?! = — 2 /?! , 



(3.22) 



The coefficients z-q, multiplying the single pole terms in the equations above, 
originate from the operators O in Eq. (2.1) and Oi in Eq. (2.2) and they are 
removed via the operator renormalization constants Zq in Eq. (2.21) and 
Zu in Eq. (2.24). 



S9,B 



X 



2\ £/2 



91 



W 



(1) 
OT,B 



m 



2\ ^ 



^ ( + I pS' ® ^i? + I P^' ® - 5 ^-i?) 



(1) 



® -Sb + + 4°^) ^ -Sb + -Sb 



(3.23) 



2\ e/2 



(1) , -(1) 



2\ s 



^ Pi? ® Pi? + J ( - 6 /3o u.^;-^ + 2 Pi? « 
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+2 P^^ ® w^l^ - 6 /3o + 2 



+2P4°)(8)w;;;-'r + w, 



qq,B ^ ^qq,B 



(3.24) 



qiq2,B ) 



-2 99 



99 



+ '^P^>^qq,B + wZ,B 



where ® denotes the convolution symbol defined by 



(3.25) 



/ ® 9{z) = / / dz2 S{z - zi Z2) f{zi) g{z2) . (3.26) 
JO Jo 

For identical quark-quark scattering we have the same formula as in Eq. 
(3.25) except that w]^^_^ ^ '"^giL.B- The expressions above follow from the 
renormalization group equations. They are constructed in such a way that 
become finite after coupling constant renormalization, operator renormaliza- 
tion and mass factorization are carried out. The splitting functions Pab{z) 
and the coefficients Wab{z) with z = m^/s also occur in the coefficient func- 
tions given below except for the NLO terms w^abi^)-> which are proportional 
to 6. They are given by 



w 



(1) 

99>H 



8V2(z) -6C{2)Vo{z) + (^--16 + 8^-8^2^ (in^(l-z) 



+—z-—z^ + 25(l-z) 
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0-4 + 2z) (ln^(l-^) -\nz\u{l- z) + ^\u^ z 



-iC(2)) + (-^ + 6-.) (ln(l-.)-iln.) + l-5 + 



-,(1) 



{l-zf 



1 1 

(ln(l-^)- ln^+-) 



3 z 

where denotes the distribution 

^ (hi'{l-z) 



1-z 



(3.27) 



(3.28) 



The difference between scalar and pseudo-scalar Higgs production is given 

by 



-8(1-^) -145(1-^) 



W 



(1) 

qq,A-R 



0, 







(3.29) 



To render the partonic cross sections finite one has first to perform coupling 
constant renormalization. This is done by replacing the bare coupling con- 
stant by the renormahzed one see Eq. (2.25). Then one has to carry out 
operator renormalization which is achieved by multiplying the expressions 
in Eqs. (3.22)-(3.25) with the operator renormalization constants in Eqs. 
(2.21), (2.24). The remaining divergences, which are of coUinear origin, are 
removed by mass factorization 



V ^ / c,d=q,q,g 



r 



db 



-1 A / "^'^ 

Z Aed,B ( -J^ 



(3.30) 



where Tca{z) denote the kernels containing the splitting functions which mul- 
tiply the coUinear divergences represented by the pole terms 1/e'^. Note that 
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in the case of pseudo-scalar Higgs production the Zq (see Eq. (2.21)) in the 
above equation has to be replaced by Z\\ in Eq. (2.24). For the four different 
subprocesses the mass factorization relations in Eq. (3.30) become equal to 



Tgg ® Tgg ^ ;2 " ^ Agg^B + 4 ^ T,, ^ ^ A 



99 ^ '- qg 



Zl W, 



O ^"^99,^ 



r„„ ®V„„® Z ^ \q,B + Tgq®Vgg® Z ^ A 



qqyy gg 

+r, 



qq ^ qg 



r„„ ® z-^ A, 



qq,B , 



Zl W, 



O ^U9^ = ^gq^^gq 



®V„„®Z ^ Agg^ + Vgq ^ V qg ® Z ^ A, 



+ 2 Vgg ® Tgg ® Z^^ A 



qq,B 



(3.31) 



(3.32) 



(3.33) 



z-^ZlW, 



O >'»'qiq2,B — 



-Egg ® Egg 



Z-lAgg,B + 2E 



gq ® 



;2-^A 

11 11 

where we have identified 

E = E- E = E - 

^99 ^99 1 ^99 ^99 1 

W^Sg,B = Wgg^B = Wgg^B = W ^g , 



z-^A 



ff9,B 



(3.34) 



^?9,B 



^?5,B 



(3.35) 



Since we need the finite expressions up to order a1 it is sufficient to expand 
the kernels Tea up to the following order in the renormalized coupling constant 



99 



5(1 - z) + a^iiJ?) S, 
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p(0) 
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(3.39) 
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5{l-z) + as{fi')S, 



p(0) 
99 



£2 V2 



.1 p{0) ^ p(0) ^ ^ p(0)\ ^ 1 p(l) 
2 99 ^ 99 ~ 99 J ~ 2£ 



(3.40) 



(PS denotes pure-singlet). After renormalization and mass factorization the 
coefficient functions have the following algebraic form 



z-^Agg^B = 5{1 - z) + asifji^) 



{5 ^-S' ® ^S' + i n'? ® n'J' -5/30 PS' + 3 A=} w ^ 
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In Appendix A and B we give the explicit expressions for the coefficient 



functions so that one can determine the coefficients P, 



(k) 



ab 



and wISJb- 



Our 



results which are expressed in the colour factors Ca, Cp and Tj agree with 
those published in [20], [19], [21] for N = 3. In the representation of the 
coefficient functions above we have put the renormalization scale fj,r equal 
to the mass factorization scale /i. If one wants to distinguish between both 
scales one can make the simple substitution 



l + ^/?oln4 



47r 



(3.45) 
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4 Total cross sections for the process 

In this section we will present total cross sections (see Eq. (2.1)) for Higgs- 
boson production in proton-proton collisions at the LHC and in proton-anti- 
proton collisions at the TEVATRON. The cross section can be written in two 
ways (for the definitions see [16], [44]) 



O'tot 



E f dy^^,iy,^^^)A^J^/pi , (4.1) 



=q,<i,9 

,2 



8{m 

where x — m^/S and is the momentum fraction luminosity defined by 

^ab{y,n')= ['-fa{u,n')fJ^,n') . (4.2) 
Jy U \U y 

However Eq. (4.1) can be also cast in the form 



with (^ab{y,i^^) =y''^ab{y,i^^), (4.3) 

where ^ab denotes the parton luminosity. For the exact cross section the ex- 
pressions in Eqs. (4.1) and (4.3) lead to the same results but if one wants to 
study the various approximations in the literature like the soft-plus- virtual 
(SV) gluon approximation then it makes a difference if Aab{z) or A.ab{z) 
will be expanded around z — 1. This difference will show up in the less 
singular terms as we will show at the end of this section. In the subsequent 
part of this section we study the dependence of the cross section on the input 
parameters like the QCD scale A, the renormalization/factorization scale jj,, 
the mass m of the Higgs boson and the input parton densities. We are also 
interested in the theoretical uncertainty of the NNLO cross section. One 
uncertainty is due to the corrections which show up beyond NNLO. They 
can be guessed from the rate of convergence represented by the fT-factor and 
the variation with respect to the scale /x. We will show that the coefficient 
functions A.ab{z) in Eq. (4.1) are dominated by the logarithmic terms of the 



25 



type ln*(l — 2;). A determination of these logarithms is of importance to 
estimate the higher order corrections. Another uncertainty is caused by the 
input parton densities. Parton density sets differ from each other in several 
aspects. They are based on fitting data from different experiments. Fur- 
thermore the shapes of the densities at the specific input scale chosen for fj, 
differ from each other. Finally these sets arc based on different theoretical 
assumptions, for instance in the treatment of heavy flavours. In some sets the 
heavy flavour is treated as a massless particle which is described by a parton 
density similar to the treatment of the light quarks. In other sets the mass of 
the heavy flavour is considered to be on the same order of magnitude as the 
other large scales. Then the heavy flavour production is described by exact 
perturbation theory at a certain order so that threshold effects are fully taken 
into account. Although the coefficient functions in the effective Lagrangian 
approach are computed exactly in NNLO this is not the case for the parton 
densities because the exact three-loop splitting functions (anomalous dimen- 
sions) arc not known yet. Until now only a finite number of moments are 
available (see [45]) which are used in [46] together with other constraints to 
approximate the splitting functions. These approximations are very reliable 
as long as y > 10~^ in Eq. (4.1). The approximated splitting functions 
were used in [47] to obtain NNLO parton density sets. One of them called 
MRST02 (see Table 1) will be used in our paper. For the approximations in 
LO and NLO we use the sets in [47] and [48] respectively. For the LO, NLO 
and NNLO plots we employ the one-, two-, and three-loop asymptotic forms 
of the running coupling constant as given in Eq. (3) of [49]. In order to 
make a comparison with other parametrizations for the parton densities we 
also choose the sets made by CTEQ [50] and GRV [51]. However these sets 
do not contain NNLO versions so that we can present the NLO results only. 
All sets are fisted in Table 1 ^ together with the corresponding QCD scale 
A„j, determined for n/ = 5. The same number of flavours is also chosen for 
the coefficient functions. Notice that the GRV sets do not contain densities 
for charm and bottom quarks. For simplicity the factorization scale is set 
equal to the renormalization scale iir- For our plots we take /x = m and 
use the MRST sets for the LO, NLO and NNLO computations (see Table 1) 
unless mentioned otherwise. Here we want to emphasize that the magnitude 

^Notice that according to the prescription of the CTEQ group [50] also the LO corrected 
quantities have to be computed with the NLO running coupHng constant. 
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MRST02 (LO, lo2002.dat) 


A^^ = 


167 MeV 




-- 0.130 


MRSTOl (NLO, alfll9.dat) 




= 239 MeV 




= 0.119 


MRST02 (NNLO, vnvalfll55.dat) 


^NNLO 


= 176 MeV 




= 0.115 


CTEQ6 (LO, ctcq61.tbl) 


^NLO , 


= 226 MeV 




= 0.118 


CTEQ6 (NLO, cteq6m.tbl) 




= 226 MeV 




= 0.118 


GRV98 (LO, grv981o.grid) 


ALO = 


132 MeV 




= 0.125 


GRV98 (NLO, grvnlm.grid) 




= 174 MeV 




= 0.114 



Table 1: Various parton density sets with their values for the QCD scale A 
and as{Mz). 

of the (Ttot is extremely sensitive to the renormalization scale because the 
effective coupling constant Gb ~ asifJr), which implies that cr^° ~ a^. The 
sensitivity to the factorization scale is much smaller. 

For the computation of the effective coupling constant Gb in Eq. (2.3) we 
choose the top quark mass rrit — 173.4 GeV/c^ and the Fermi constant Gp — 
1.16639 GcV~^ = 4541.68 pb. In this paper we will only study scalar Higgs 
production and omit pseudo-scalar Higgs production. The cross sections of 
the latter are about 9/4 cot^ (3 larger than those for the standard model Higgs 
so that all our conclusions also apply to pseudo-scalar Higgs production. We 
give results for both proton-proton collisions at the centre-of-mass energy 
\/S = 14 TeV and proton-anti-proton collisions at y/S = 2 TeV. 

In Fig. 1 we have plotted the contributions coming from the various sub- 
processes up to NLO. From this figure we infer that (Xtot is completely dom- 
inated by the gg reaction whereas the other contributions are down by two 
{gq + 51^-subprocess) or even by three (g'^-subprocess) orders of magnitude. 
This picture is different from the behaviour of the transverse momentum p^- 
and rapidity ^/-distributions in [14] where the gq + (jfg-reaction amounted to 
about 1/3 of that for the gfgf- reaction when px > 30 GeV/c. This is because 
(Ttot receives all its contribution from the small pr-region {px < 30 GeV/c) 
where the gg subprocess overwhelms all other reactions. This picture is not 
changed when wc study the cross section in NNLO (see Fig. 2) where new 
subprocesses contribute given by the qq, qq and qq reactions. Notice that 
the incoming (anti-) quarks can be identical and non-identical. The relative 
order of magnitude is the same as in NLO and the new reactions are down 
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by three orders of magnitude. In Fig. 3 we have plotted o"tot in LO, NLO 
and NNLO as a function of m. The cross section decreases when m increases 
due to a reduction in the available phase space. In NNLO the curve shows a 
steeper decrease than in NLO and LO. Furthermore the growth of the cross 
section slows down while going to higher orders. To show more clearly how 
the cross sections change with respect to m we have put them in Tabic 2. 
We have also shown them for the TEVATRON (proton-anti-proton collisions 
at y/S — 2 TeV) in Table 3 where they are appreciably smaller than in the 
case of the LHC. 

The NNLO cross sections in Tables 2, 3 have some theoretical errors which 
can be estimated in various ways. First one can study their variation with 
respect to the scale /x. This can be achieved by computing the ratio 




where = m and // is varied in the range 0.1 < /i/jio < 10. In Fig. 4 a plot 
of N is shown for m = 100 GeV/c^. Here one observes a clear improvement 
while going from LO to NNLO. In particular the curve for NNLO is flatter 
than that for NLO. The improvement becomes even better when we choose 
a larger Higgs mass see e.g. Fig. 5 {m — 200 GeV/c^). However for still 
larger masses there is hardly any improvement anymore (see Fig. 6 where 
m = 300 GcV/c^). This is in contrast to the NNLO corrected cross section 
for vector boson production in [26] which is insensitive to the choice of scale 
chosen in the range above. 

A second way to study the reliability of our prediction is to study the 
rate of convergence of the perturbation series, which is represented by the 
X-factor. We choose the following definitions 



^NLO _.NNLO 
Tv-NLO "tot Tv-NNLO "tot /a c-\ 

"tot "tot 



Notice that in the above expression the definition of i^^'^^'^o jiff^j^g from the 
usual one which is given by X^nlo _ ^^nnlo^'^lo hsive chosen this 

definition because the rate of convergence can be shown in a better way. 
In Fig. 7 one observes that both X-factors vary slowly as m increases. 
Moreover there is considerable improvement in the rate of convergence if 
one goes from NLO to NNLO. At m = 100 GeV/c^ we have K^^^ ~ 1.74 
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mass 


LO 


NLO 


NNLO 


100 


30.35 


52.75 


60.84 


110 


25.53 


44.75 


51.80 


120 


21.77 


38.43 


44.62 


130 


18.77 


33.37 


38.85 


140 


16.36 


29.27 


34.17 


150 


14.38 


25.88 


30.29 


160 


12.74 


23.06 


27.07 


170 


11.37 


20.69 


24.36 


180 


10.23 


18.69 


22.09 


190 


9.254 


17.00 


20.10 


200 


8.425 


15.53 


18.43 


210 


7.714 


14.28 


16.97 


220 


7.102 


13.20 


15.71 


230 


6.573 


12.26 


14.62 


240 


6.115 


11.45 


13.69 


250 


5.718 


10.74 


12.86 


260 


5.375 


10.13 


12.16 


270 


5.080 


9.604 


11.55 


280 


4.827 


9.152 


11.03 


290 


4.615 


8.772 


10.59 


300 


4.441 


8.468 


10.24 



Tabic 2: Total cross sections in pb for Higgs masses between 100 GcV and 
300 GcV at the LHC. The LO. NLO and NNLO results are generated with 
the MRST parton densities listed in Table 1. 
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mass 


LO 


NLO 


NNLO 


100 


0.6116 


1.351 


1.781 


110 


0.4608 


1.025 


1.363 


120 


0.3530 


0.791 


1.060 


130 


0.2745 


0.619 


0.835 


140 


0.2161 


0.491 


0.666 


150 


0.1721 


0.393 


0.538 


160 


0.1385 


0.318 


0.438 


170 


0.1124 


0.259 


0.359 


180 


0.0920 


0.213 


0.293 


190 


0.0758 


0.177 


0.248 


200 


0.0629 


0.147 


0.207 


210 


0.0526 


0.124 


0.176 


220 


0.0442 


0.105 


0.149 


230 


0.0374 


0.089 


0.127 


240 


0.0318 


0.076 


0.109 


250 


0.0270 


0.065 


0.095 


260 


0.0234 


0.056 


0.082 


270 


0.0202 


0.049 


0.072 


280 


0.0176 


0.043 


0.063 


290 


0.0154 


0.038 


0.055 


300 


0.0136 


0.033 


0.049 



Tabic 3: Total cross sections in pb for Higgs masses between 100 GeV and 
300 GeV at the TEVATRON. The LO, NLO and NNLO results are generated 
with the MRST parton densities listed in Table 1. 
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whereas i^^^^o ^ 1.15. gtiH the corrections for Higgs boson production are 
larger than those obtained from Z-boson production at the LHC where one 
gets X^LO ~ 1.22 and ^nnlo ^ q.OS (see [26]). In Fig. 8 we compare 
the X-factor in NLO obtained from the MRST with the results computed 
by using the parton density sets given by the CTEQ [50] and GRV [51] 
collaborations (see Tabic 1). From this figure wc infer that the X-factor 
computed using GRV98 and CTEQ6 is larger than the one obtained from 
MRSTOl(NLO) at m = 100 GeV/c^ However at large m there is a cross 
over and we observe /^cteq ^ ^mrst ^ j^grv98 ^ rjj-^g large result shown 

by the CTEQ6-sct is due to the behaviour of their gluon density in LO. In 
NLO this density behaves in the same way as those in the MRST and GRV 
sets. We also investigated older sets provided by the GTEQ-coUaboration, 
namely CTEQ5 and CTEQ6. It turns out that all these LO CTEQ sets lead 
to the same large value for K^^'^ as shown by Fig. 8. A study of K^^*^, 
using the most recent parton densities, was also made in [44] but now for the 
TEVATRON. 

After the K-factor we now make a comparison between cxtot computed 
from the various sets which leads to the third uncertainty in the theoretical 
prediction. For that purpose we plot the ratios 

^CTEQ GRV 

pCTEQ _ ^tot pGRV _ "tot (a ^\ 

~ ^MRST ' ~ ^MRST ' \^-^) 

^tot CTtot 

which are shown in NLO in Fig. 9. In this figure we see that the CTEQ6 set 
leads to the same cross section as the one obtained by MRSTOl. The cross 
section given by the GRV98 is 20% above the MRSTOl result and the ratio 
decreases with increasing m to 5% above at m = 300 GeV/c^. The reason 
that the GTEQ6 result is closer to the MRST sets can be attributed to the 
fact that their parton density sets are more recent and are constructed to fit 
the same experiments. Notice that a study of the ratios in Eq. (4.6) was also 
made for the TEVATRON in [44] using the most recent parton densities. 

The factorization scale dependence in Figs. 4-6 and the /^-factor in Fig. 
7 of the NNLO cross section can be used to give an error estimate of the 
latter for the LHC. Since ((j^^^o _ ^nlo^/^nlo ^ ^nnlo _ ^ ^^^^^ 
^NNLO jg gq^j^i ;l5%, 19% and 21% for m = 100 GeV/c^ m = 200 GeV/c^ 
and m = 300 GeV/c^ respectively (see Fig. 7). This is corroborated by the 
scale variation if we choose the range 0.25 < fi/fiQ < 4.0. For fi/ fiQ = 0.25 the 
errors above become 18%, 15% and 14% respectively. In the case n/no = 0.25 
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they are given by 16% for all m given above. Notice that the error due the 
i^-factor increases when m gets larger whereas the opposite happens if the 
error estimate is derived from the mass factorization/renormalisation scale. 
On top of this we have to add the error due to the chosen parton density. For 
this we choose the NLO plots in Fig. 9 where the cross section predicted by 
CTEQ6 hardly differs from the one given by MRSTOl but the GRV98 plot 
deviates from the latter by 20% for small m to 5% for large m. 

An important feature of the Higgs cross section, which is very often em- 
phasised in the literature, is that the integral in Eq. (4.1) is dominated by 
y X = VP? j S. Since x is small the gg part of o"tot dominates the contribu- 
tions coming from the other partonic subprocesses because of the steep rise 
of the gluon-gluon luminosity in Eq. (4.2) at small y. The importance of the 
small y-region is revealed when we impose an upper cut on the integral in 
Eq. (4.1). Hence we compute 



'-'^tot('^max) — 



dy^ab{y,l^ ) 

a,b=q,g i=l ^ 



xA 



and plot the ratio 



ab 



X m 



-^(•^max) 



''"tot (-^max) 
O'tot(l) 



(7tot(l) 



EXACT 
'tot 



(4.7) 



(4.8) 



which is shown for NLO and NNLO in Fig. 10a for the choice Xmax = 
5 X. The figure reveals that more than 95 % of the cross section comes 
from the integration region x < y < 5 x where x < 5.1 x 10~^. When m 
gets larger R{x^ax) will increase because the available phase space for Higgs 
boson production decreases. The latter also happens when the centre- of- mass 
energy decreases like in the case for the TEVATRON where = 2 TeV 
(Here x < 2.5 x 10~^) see Fig. 10b. (Note Figs. 10a and 10b have different 
scales.) In the latter case more than 99 % of the cross section receives its 
contribution from a; < y < 5 a; which becomes 100 % when m — 300 GeV/c^. 

Finally we study the validity of the soft-plus- virtual gluon approximation 
including subleading terms. The result of this approximation depends on the 
definitions for the total cross section given in Eq. (4.1) and Eq. (4.3) as has 
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been pointed out in [44]. If we follow the definition in [8], [17], where one 
choses the expression in Eq. (4.1), the sof-plus- virtual coefficient function 
becomes 




i=l L j=0 



2i-l 



H ai,jVj + bi5{l- z) 



.(4.9) 



where the distributions T>i{z) are defined in Eq. (3.28). Furthermore the 
logarithmic terms ln-'(l — z) including the constants {j = 0) are important. 
The latter are denoted by 
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Notice that the coefficients aij,b.i,Cij also contain terms in Inm^///^. The 
soft-plus- virtual gluon contributions only emerge from the (jf^f-subprocess and 
they can be found in Eqs. (A. 2) and (A. 6). This subprocess also leads to 
the logarithms of the type W(1 — z). The latter also show up in the 
gq + gq-channel but they are absent in the qq and qq subprocesses provided 
one neglects the terms (1 — z)'' for k > 0. The coefficients q j can be found 
in Eqs. (A.26)-(A.32) (see also [18]) (scalar Higgs) and Eqs. (B.21), (B.22) 
(pseudo-scalar Higgs). If we use the definition of the cross section in Eq. 
(4.3) advocated in [16], [44] the expansion of z~^ A^^ will not differ from the 
one shown in Eq. (4.9) but A^^ will change due to contribution of subleading 
logarithms and it becomes equal to 
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E<(/^') E + c«) - ^) ■ (4-11) 

\ A* / i=l j=0 

Because of the different luminosities ^{z) and ^{z) in Eqs. (4.1) and (4.3) 
respectively the value of the soft-plus-virtual cross section af^ will change. 
This will be only partially compensated by al'^^ coming from Eqs. (4.10) and 
(4.11). In fact the two results for crf^'" coming from definitions (4.1) and 
(4.3) will differ by terms in (1 — z)'^ ln-'(l — z) for k > 1 convoluted by one 
of the luminosities. In order to study the effect of the two definitions on the 
approximation we will plot 

SV SVL 
Ctot CTtot 
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The same study is made in [8], [16], [17] and [44]. However in these papers 
the complete result for the NNLO coefficient function was not known yet. 
Therefore one could only study the leading logarithmic terms in Eqs. (4.10), 
(4.11) which are of colhncar origin. In Fig. 11a and Fig. lib we have shown 
the ratios above for the LHC in NLO and NNLO respectively. The figures 
reveal that af^^ using Eq. (4.3) gives a much better approximation than Eq. 
(4.1). This also holds if we include the logarithms in Eqs. (4.10) and (4.11) 
denoted by crf^^. The latter overestimates the exact cross section a httle bit. 
Further we observe that both SV approximations become worse in NNLO 
with respect to NLO. In Figs. 12a,b we have performed the same analysis 
but now for the TEVATRON. The features are the same as in Figs. lla,b 
for LHC except that for the TEVATRON the SV approximations derived 
from Eq. (4.1) and Eq. (4.3) become better. This is no surprise because for 
the TEVATRON the Higgs mass is much closer to the boundary of phase 
space. This also explains why at larger masses the SV approximation im- 
proves. From this study we infer that the soft-plus-virtual approximation 
works much better for Eq. (4.3) than for Eq. (4.1). The reason is that the 
cross section is dominated by the small x- region (see Eq. (4.7) and Fig. 10). 
By choosing the gluon luminosity in Eq. (4.3) this region will become even 
more important than the gluon luminosity chosen in Eq. (4.2) which is used 
in Eq. (4.7). For NLO this was already shown in [44] but it is now also con- 
firmed in NNLO. Further this approximation becomes worse in NNLO which 
holds for both definitions. In spite of this drawback one should still choose 
Eq. (4.3) to perform the soft-plus- virtual gluon resummation (see [52]) in 
order to get a better estimate of the cross section. We also noticed that the 
contribution of logarithmic terms in Eqs. (4.10), Eqs. (4.11) is substantial. 
The determination of the coefficient of the leading term, which is of coUinear 
origin, is done in [8], [16]. A systematic approach to determine can be 
found in [53]. 

Summarising our results we have recalculated the NNLO corrections to 
the total cross section for (pseudo-) scalar Higgs production using a different 
method than those presented in [20] and [21]. In particular our approach dif- 
fers from the one given in [20] since we directly evaluated the multi-particle 
phase space integrals without fitting the coefficient functions to an expansion 
in the terms {1 — z)'' \n\l — z) . Further we presented the radiative parts of 
the coefficient functions for general colour factors and checked that for = 3 
our answers are in agreement with the literature. It turns out that the total 
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cross section is almost completely determined by the gg-subprocess which is 
in contrast to the differential cross section where also the gq + ^fg-channel 
gives a considerable contribution. We studied the X-f actors and the depen- 
dence of the cross sections on the chosen mass factorization/renormalization 
scale and the adopted parton density set from which one can infer an er- 
ror estimate of the cross section. If we exclude the GRV98 parton density 
the error can be mainly attributed to the missing higher order contributions 
to the coefficient functions which we estimate to lie between 14% and 21%. 
Finally we studied the validity of the soft-plus- virtual gluon approximation. 
Depending on the definition of the total cross section (Eq. (4.1) versus Eq. 
(4.3)) this approximation is excellent in NLO but it becomes less good in 
NNLO but Eq. (4.3) is the best way to resum the soft-plus-virtual correc- 
tions. 
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Appendix A 



In this Appendix we present the coefficient functions for scalar Higgs boson 
production. In the case of the gg subprocess one can split the corresponding 
coefficient function into 



(A.l) 



where the superscripts S + V and H denote the soft-plus- virtual and hard 
gluon parts respectively. The former contains the distributions 5{1 — z) and 
'Di{z) where the latter is defined in Eq.(3.28). We find 
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Notice that the colour decomposition of the S{1 — z) term into UfTfCA and 
UfTfCp is arbitrary due to our ignorance of the same parts appearing in the 
two- loop virtual corrections computed in [15]. They are only correct when 
N = 3. 
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As discussed below Eq. (4.10) the corrections are dominated by the terms 
Cij ln-'(l — z) in the coefficient functions. They are obtained by taking the 
limit ^ — > 1 in the expressions above. Most of them vanish. Those which 
survive are given by 
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The nonvanishing coefficients in Eq. (4.9) can be read off from tlie equa- 
tions above. 
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Appendix B 



In this Appendix we present the coefficient functions for pseudo-scalar Higgs 

boson production. To shorten the expressions we give the difference between 
the coefficient functions originating from pseudo-scalar and scalar production 
only, namely, 
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In the NNLO corrections below we indicate by O12 the contribution coming 
from the interferences between the operators Oi and O2 in Eqs. (3.9), Eqs. 
(3.10). In the coefficient functions one has to put O12 = 1, 
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For the colour decomposition of the n/-part of the expression above see the 
remark below Eq. (A. 6). Also 
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The nonvanishing coefficients in Eq. (4.9) appearing in the difference 
between the pseudo-scalar and scalar coefficient functions are determined 
from the expressions above by taking the limit 2; ^ 1. Those which survive 
in this limit are given by 
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Figure Captions 



Fig. 1 The total cross section a^ot plotted as a function of the Higgs mass at 
= 14 TcV with jj, = m. The NLO plots are presented for the total 
(solid line) and the subprocesses gg (long-dashed line), abs{gq+ gq) 
(dot-dashed line) and 100 x {qq) (short-dashed line). Also shown is the 
LO result (dotted hne). 

Fig. 2 The total cross section cxtot plotted as a function of the Higgs mass at 
^fS = U TeV with jj, = m. The NNLO plots are presented for the total 
(solid line) and the subprocesses gg (long-dashed line), Sihs{gq + gq) 
(dot-dashed hne), 100 x [qq) (dotted line) and 100 x [qq -\- qq) (short- 
dashed hne). 

Fig. 3 The total cross section dtot with all channels included plotted as a 
function of the Higgs mass at \fS = U TeV with II — m; NNLO (solid 
line), NLO (dashed line) and LO (dotted line). 

Fig. 4 The quantity N{h/hq) (see Eq. (4.3) at = 14 TeV plotted in the 
range 0.1 < fi/fiQ < 10 with fiQ = m and m = 100 GeV/c^. The results 
are shown for LO (dotted line), NLO (dashed hne) and NNLO (solid 
line) . 

Fig. 5 Same as Fig. 4 but now for m = 200 GeV/c^. 
Fig. 6 Same as Fig. 4 but now for m = 300 GeV/c^. 

Fig. 7 The iT-factors in NLO and NNLO at \/S = U TeV as a function 
of the Higgs mass using the MRST-sets; /sT^nlo ^^^^^^ ^^^^^^ j^nlo 
(dot-dashed line). 

Fig. 8 The A'-factors in NLO at = 14 TeV as a function of the Higgs 
mass using the following parton density sets MRSTOl (solid line), 
GRV98 (dashed hne) and CTEQ6 (dotted line). 

Fig. 9 The ratios R = atot/^t^t^^^ (Eq.(4.5)) in NLO at = 14 TeV 
and /i = m as a function of the Higgs mass; R^'^^ (sohd hne), T^^teq 
(dotted line). 
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Fig. 10a The ratio -R(.'riiiax) (sec Eq. (4.7)) for proton-proton collisions 
(LHC) where Xmax = 5 a; with x = w?/S. The CM energy and scale 
are given by — 14 TeV and ji — m respectively; NNLO (sohd line), 
NLO (dashed hne). 

Fig. 10b Same as in Fig. 10a but for proton-anti-proton collisions (TEVA- 
TRON) at = 2 TeV and // = m. 

Fig. 11a The ratios in Eq. (4.12) in NLO for proton- proton collisions (LHC) 
as a function of the Higgs mass at \/S — 14 TeV and ji — m; it!^^(Eq. 
(4.1)) (dot-dashed line), i?svL(Eq. (4.1)) (dashed line), i?sv(Eq. (4.3)) 
(dotted line), i?svL(Eq. (4.3)) (solid line). 

Fig. lib The ratios in Eq. (4.12) in NNLO for proton-proton collisions 
(LHC) as a function of the Higgs mass at = 14 TeV and n — 
m- i?sv(Eq. (4.1)) (dot-dashed line), R^"^^{Y.(i. (4.1)) (dashed line), 
i?sv(Eq. (4.3)) (dotted line), i^^VL^gq ^4 ^^^^id line). 

Fig. 12a Same as in Fig. 11a but for proton-anti-proton coUisions (TEVA- 
TRON) at = 2 TeV and = m. 

Fig. 12b Same as in Fig. lib but now for proton-anti- proton coUisions 
(TEVATRON) at = 2 TeV and = m. 
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